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Transverse limits on the uni-directional pulse propagation approximation 
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I calculate the limitations on the widely-used forward-only (uni-directional) propagation assump- 
tion by considering the effects of transverse effects (e.g. diffraction). The starting point is the scalar 
second order wave equation, and simple predictions are made which aim to clarify the forward- 
backward coupling limits on diffraction strength. The result is unsurprising, being based on the 
ratio of transverse and total wave vectors, but the intent is to present a derivation directly compa- 
rable to a recently published nonlinearity constrained limits on the uni-directional approximation 
i- 



I. INTRODUCTION 

Most approaches to optical pulse propagation rely on 
an approximation where the fields only propagate for- 
wards. Even the recently derived extensions of typical 
propagation methods used in nonlinear optics (e.g. [3,3) 
assume a complete decoupling between oppositely propa- 
gating fields to optimize the calculation. Moreover, those 
based directly on Maxwell's equations (e^. [1, H, Q) or 
the second order wave equation (e.g. 0,13 Q), are often 
simplified to work in the forward-only limit, where back- 
ward propagating fields are set to zero. This is despite 
directional decompositions of Maxwell's equations (e.g. 

indicating that any effect not allowed for by that 
decomposition couples the forward and backward waves 
together - and even creates a backward field if one is 
not present. Usually we assume that a forward wave will 
not generate a significant backward wave via the diffrac- 
tion because we are in the paraxial limit, and even then 
any generated backward component is very poorly phase 
matched^ . 

I compare predictions for propagation wave vector 
from uni- and bi-directional theories. Such a compar- 
ison has been done for the more important case of 
nonlinearity-induced forward-backward coupling 
used analytical expressions for carrier shocking 
in concert with simulations to examine the effects of a 
uni-directional approximation. Here my intention is sim- 
ply to provide a complementary calculation to clarify the 
effects of forward-backward coupling induced by trans- 
verse (diffraction) effects. The result, while unremark- 
able, does clarify the bounds on the validity of propaga- 
tion models using a uni-directional approximation. 

Since linear dispersion and finite nonlinear response 
times will typically diminish any generation of a back- 
ward wave, it is clear that any model which can be as- 
sumed uni-directional on the basis of this paper will be 
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^ If the forward field has a wave vector fco evolving as exp(+4fcoz), 
the generated backward component will evolve as exp(— ifcoz). 
This gives a very rapid relative oscillation exp(— 2jfco2), which 
will quickly average to zero. 



more so in practice. These results do not tell us whether 
the uni-directional approximation would be more or less 
robust for situations requiring vector fields (e.g. Ts]), or 
for nonlinear effects such as self- focusing ^14i] . or nonlin- 
ear diffraction [l^ but they at least establish a point of 
reference. 



II. BASIC THEORY 

Most optical pulse problems consider a uniform and 
source free dielectric medium. In such cases a good start- 
ing point is the second order wave equation, which results 
from the substitution of the W x H Maxwell's equation 
into the V x E one in the source-free case (see e.g. [l6|). 
Further, assuming linearly polarized pulses, we can use a 



scalar form. Defining V"^ = + dy - 
we can write the wave equation as 



9^ and da = d/do 



47r . 



E{t) ^ —d PTiE{t),t). 



(1) 



Here I have suppressed the spatial coordinates for no- 
tational simplicity; in fact we have E{t) = E{t, r) and 
the total polarization P{E{t),t) = PT{E{t,f),t,f); also 
r = {x,y,z). Here we will consider only isotropic linear 
media, which enables us to replace Pt with a refractive 
index; more complicated polarization behaviour is cov- 
ered elsewhere [TtI. [l8l|. Thus, in the frequency domain, 
we can write 



V 



2 n'^{uj)u)^'^ 



E{t) = 0. 



(2) 



However, in most descriptions of pulse propagation we 
will want to chose a specific propagation direction (e.g. 
along the z-axis), and then denote the orthogonal compo- 
nents (i.e. along x and y) as transverse behaviour. Thus, 
splitting the operator into its propagation (z) and 
transverse {x, y) parts, we can rewrite the wave equation 
similarly: 



[dl + K^] E{u) = -VlE{u) 
[dl + K^] E{cj) = +klE{oj), 



(3) 
(4) 
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where the total wavevector is given by = r?uj^ 
the transverse component is fc_L. If we want to describe 
diffraction, then we can give the field some suitable beam 
profile E{x^ y), an even simpler case is that of off-axis one 
dimensional propagation, which requires merely a fixed 
value of k±. 

I now factorize the wave equation, a process which, 
while used in optics for some time j7(| has only recently 
been used to its full potential 0, U 0] ■ Factorization 
takes its name from the fact that the LHS of eqn. ([3]) is a 
simple difference of squares which might be factorizcd, in- 
deed this is what was done in a somewhat ad hoc fashion 
by Blow and Wood in 1989 Since the factors are just 
dz T iK, we can see that each (by itself) would generate 
a forward directed wave equation, and the other a back- 
ward one. Without going into detail (although see the ap- 
pendix) , a rigorous factorization procedure [1, [13] allows 
us to define a pair of counter-propagating Greens func- 
tions, and so divide the second order wave equation into 
a pair of coupled counter-propagating first order ones. 

Counter-propagating wave equations suggest counter 
propagating fields, so I split the electric field up accord- 
ingly into forward {E'^) and backward {E~) parts, with 
E = E~^ + E~ . The coupled first order wave equations 
are 



so that 



dzE^ 



±iKE^ 



2K 



(5) 



The RHS now falls into two parts, which I term the un- 
derlying and residual parts [13]. First, there is the iKE^ 
term that, by itself, would describe plane- wave like prop- 
agation. Second, the remaining part (here proportional 
to k\) which can be called "residual" terms. These resid- 
ual contributions, here containing the transverse effects, 
account for the discrepancy between the true propagation 
and the underlying propagation. Although here the resid- 
ual component will be only a weak perturbation in e.g. 
the paraxial limit, the theory presented here is valid for 
any strength. Although my preference would be to use a 
directional fields approach [6[ rather than the factoriza- 
tion one used here, it is difficult to describe transverse 
effects satisfactorally. 

Note that the work of Weston examines this kind of 
wave-splitting with more mathematical rigour (see e.g. 
[20|). although without consideration of residual terms, 
and (at least initially) in the context of reflections and 
scattering. This theory was based on that from the ear- 
lier work of Beezley and Krueger ^21] who applied wave- 
splitting concepts to optics. 



A. Bi-directional (exact) case 

The scalar second order wave equation given above in 
eqn.(IT|) trivially provides a total wavector for any given 
direction of propagation. This is simply a sum of squares 
of the parallel and transverse contributions fc|| and fcx, 



= k\ 



(6) 



Note the reversed signs between the RHS terms of 
eqn. ([5]), the transverse part retards the propagation 
given by the underlying part; the net forward-direction 
wavevector for some fc^ is therefore less than the total 
wave vector - just as would be expected. 



B. Uni-directional approximation 

Now I make the uni-directional assumption and set 
E~ = in eqn.©, so that we get a wave equation with 
underlying (oc K) and residual (oc k']_/2K) components, 
i.e. 



dzE^ = ±iKE^ 



2K 



+ 



(7) 



Note that the diffaction term here is identical to that ob- 
tained by applying the standard paraxial approximation 
to propagation in a linear dispersive medium^. 

Alternatively, since there is no E~ field to complicate 
matters, I might rewrite this wave equation using a new 
uni-directional wave vector Ku to define only an under- 
lying propagation, i.e. 



dzE^ = ±iKuE^ 



where is 



K,, = K 



2K 



= K 



1 - 



1 fci 

2 if 2 



(8) 



(9) 



As expected, Ku is not equivalent to the true wave vec- 
tor K; indeed we expect it to be an approximation to 
/c|| , which specifies on-axis spatial variation of the exact 
propagation. 



III. FORWARD-BACKWARD COUPLING 

In the above, we saw that the bi-directional and uni- 
directional models gave different propagation wave vec- 
tors. However, note that when k'^/K^ <C 1, terms of or- 
der {k±/K)'^ or higher are negligible. We can rearrange 
and then approximate eqn.® in that limit so that 



ku = K 



1 - 



if2 



K 



2^2 



(10) 
(11) 
(12) 
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Essentially what the condition k^/K"^ <^ 1 means is 
that propagation effects transverse to the chosen propa- 
gation direction must occur on a scale much larger than 
one wavelength, or else a uni-directional approximation 
will fail. 

It is important to note that the existence of significant 
forward-backward coupling does not always demand the 
presence or generation of a freely propagating backward 
wave. It is possible for the backward wave (i.e. E~) to 
be dragged along by the forward one, as seen for several 
example in the directional fields formalism of Kinsler et 
al. [y]. Nevertheless, although it such a situation might 
give an answer correct to within a suitable scaling, in such 
a case the uni-directional approximation is not strictly 
valid. 



tion. This was done by a simple comparison of wave vec- 
tors obtained from electromagnetic scalar wave equations 
allowing for all three spatial dimensions; using both an 
exact (and hence bi-directional) model, and an approx- 
imate uni-directional model. These results are done in 
the same style as, and are intended to complement ex- 
isting limits placed on nonlinear effects fH; they are not 
intended to startle the reader with their novelty. 

I have shown that the condition k±_ /K"^ <^ 1 must hold 
for the uni-directional approximation to be true; even 
when no backward field is initially present. Unsurpris- 
ingly this is comparable to the condition for the widely 
used "paraxial" limit. Note, however, that the use of the 
paraxial approximation is rarely accompanied by a dis- 
cussion of potential generation of backward propagating 
waves. 



IV. CONCLUSION 

I have demonstrated one of the fundamental limits on 
the widely used uni-directional propagation approxima- 
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Appendix: Factorizing 

Here is a quick derivation of the factorization process; 
the z-derivative has been converted to «fc, 0^ = n^uP' j(? ^ 
and the unspecified residual term is denoted Q. 



E = 



1 



k^ - (3- 
-1 



■Q 



1 



{k-P) (fc + /3) 



1 



1 



k + P k-(3 



Q. 



(13) 
(14) 

(15) 



Now {k — (3)~^ is a forward-like propagator for the field, 
and (k + /3)^^ a backward- like propagator. Hence write 
E = E^ + E^ , and split the two sides up 



E+ +E- 
E^ 

[kTH] 

ikE^ 



-1 



1 



_k + (3 
±1 1 „ 



k-(3 



Q 



2f3 kTf3 
1 1 



= ±- 



'2(3kTP 



Q 



±^PE- ± -Q 



(16) 
(17) 
(18) 
(19) 



and reverting to z derivatives gives us the final form 



d,E^ 



±ipE^ ±—Q 
2p 



(20) 
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